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Abstract The dimer-atom-atom recombination process in the system of four iden- 
tical bosons with resonant interactions is studied. The description uses the exact 
Alt, Grassberger and Sandhas equations for the four-particle transition opera- 
tors that are solved in the momentum-space framework. The dimer-dimer and 
atom-trimer channel contributions to the ultracold dimer-atom-atom recombina- 
tion rate are calculated. The dimer-atom-atom recombination rate greatly exceeds 
the three-atom recombination rate. 
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1 Introduction 

Few-particle systems with large two-particle scattering length a possess universal 
properties that are independent of the short-range interaction details. The three- 
particle system was investigated theoretically by V. Efimov more than 40 years 
ago [T] but only in the last decade the cold-atom physics experiments [2J confirmed 
his prediction for the existence of zero orbital angular momentum weakly bound 
three-particle states with asymptotic discrete scaling symmetry. This boosted the 
interest also in the universal systems with four or even more particles. In contrast 
to the three-body system where semi-analytical results have been obtained (see 
Ref. [3| for a review), most of the theoretical studies of the four-body systems are 
numerical. In addition to the numerous bound-state calculations, e.g., [ll[5l[6l[71[8l 
l9l ll0llll| . also the collision processes have been investigated in the framework of 
hyperspherical harmonics (HH) [6], coordinate-space Faddeev-Yakubovsky (FY) 
equations [12|[8] or momentum-space Alt, Grassberger and Sandhas (AGS) equa- 
tions [131114] . The studies include the elastic and inelastic atom-trimer (8l ll4|[T5] 
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Fig. 1 (Color online) Schematic representation of the two- and three-cluster thresholds in 
the four-boson system as functions of the two-boson scattering length. The intersection of the 
dimer-atom-atom (solid line) and the nth atom-trimer (dashed line) thresholds defines aj^. 



and dimer-dimer scattering [161117] as well as the four-atom recombination [18ll6l 
1191120] , Depending on the collision energy and the two-particle scattering length 
(that can be manipulated by the external magnetic field in the cold-atom ex- 
periments) all these reactions may get resonantly enhanced due to the presence of 
unstable tetramer states. Such a resonant behaviour in the dimer-dimer relaxation 
and four-atom recombination, roughly consistent with theoretical predictions, was 
observed recently in experiments with ultracold atoms performed near the univer- 
sal regime [2n[2il[23l[23] . 

The four-particle reaction, that, to the best of our knowledge, has not been 
studied so far, is the three-cluster recombination, i.e., an inelastic collision of a 
dimer and two atoms leading to the two-cluster final state either with two dimers 
or with one atom and a trimer. Due to time reversal symmetry, the amplitude for 
this process is equal to the amplitude for the three-cluster breakup of the initial 
two-cluster state. The latter is more directly related to the transition operators for 
the two-cluster reactions ^14ul7u20) and is derived in the present work using the 
formalism of the AGS four-particle scattering equations. Their numerical solution 
is performed in the momentum-space framework. 

We will study the dimer-atom-atom recombination in the system of four iden- 
tical bosons with resonant two-boson interactions. In the potential model taken 
over from Refs. [14111 TIED] there is only one weakly bound dimer but many trimers 
since the system exhibits an Efimov effect [ll[3]- In our nomenclature we label 
the Efimov trimers by one integer number n, starting with n = for the ground 
state. The dependence of various scattering thresholds on the two-boson scattering 
length a is schematically shown in Fig. [T] The special value of a corresponding to 
the intersection of the n-th atom-trimer and dimer-atom-atom threshold is denoted 
by ai. 

In Sec. [2] we derive the three-cluster breakup amplitude in terms of the AGS 
transition operators and discuss some technical aspects of the calculations. In 
Sec. [3] we present results for the dimer-atom-atom recombination. We summarize 
in Sec. H 
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2 Four-particle scattering 

We consider the relative motion in the four-particle system interacting via short- 
range pairwise potentials vj where j labels one of the six pairs. The full Hamilto- 
nian is 

6 

H = Ho + J2^i (1) 

with Hq being the kinetic energy operator. We do not calculate the four-particle 
wave function directly by solving the Schrodinger equation or the Faddeev-Yakubovsky 
equations [12| . We prefer the integral form of the scattering equations as formulated 
by Alt, Grassberger and Sandhas 13 for the four-particle transition operators W^p. 
These equations are 18-component matrix equations 13,25 , i.e., 

i^ap = {Goti Go)^^ Sap Sji + ^ Sa^Uil'GQ tk GoUjp, (2) 

~fk 

where the components are distinguished by the chains of partitions, i.e., by the 
two-cluster partition and by the three-cluster partition. The two-cluster partitions, 
denoted by Greek letters, are either of 3 -|- 1 or 2 + 2 type. The three-cluster 
partitions, denoted by Latin letters, are of 2 + 1 -|- 1 type and therefore are fully 
characterized by the pair of particles in the composite cluster. Obviously, the pairs 
i, j and k must be internal to the respective two-cluster partitions, i.e., i G p, 
j, fc C 7 and j C cr. Furthermore, Sap = 1 — S^p and 

Go = {E + iO ~ Hor^ (3) 

is the free resolvent at the system energy E. The two-particle interactions for each 
pair are summed up to all orders to form the respective two-particle transition 
matrices 

tj = Vj + VjGotj, (4) 

while all the interactions within each two-cluster subsystem lead to the transition 
operators 

Uiif'' = Gq^ +'ySjitiGoUj'^ . (5) 

i 

The wave function |!f'p,n) describing the scattering process initiated by the 
asymptotic two-cluster channel state |<Pp,n) with the relative two-cluster momen- 
tum Pp,n is determined by the AGS transition operators as 

\1'p,n) = \'l'p,n) + J2^0tl Go Ul^^Go tk GoU^]>Wp,n)- (6) 

The energy parameter in the operators of Eq. ([6} is _E = ep,„ + p'^.n/'^P'P with ep^n 
being the energy of the nth bound state in the partition p and /ip being the reduced 
mass. The asymptotic channel state \$p,n) is a product of the nth bound state wave 
function in the partition p and the plane wave with momentum Pp.n between the 
two clusters. \'Pp,n) = \4>p,n) is decomposed into its Faddeev components that 
are calculated from the Faddeev equation 

K.n) = GoY^S^jWp.n) (7) 
j 
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and normalized such that (<?o-,ra'l^p,n) = ^crp ^n'n ^{Pa.n' ~ Pp,n)- 

The scattering amplitudes for the reactions with initial and final two-cluster 
states are given by the on-shell matrix elements of the AGS operators J^^j {(f>-'^ ^, \Ua''p\(f)p^n) 
as explained in Refs. [131125] . The amplitude for the three-cluster breakup of the 
initial two-cluster state can be obtained from Ua\\4>^p.n) as well. We derive it from 
the general relation 

6 

{^j\Tjp\'^p,n) = {1>j\Y^5,,V^<I'p^n). (8) 

i=l 

The asymptotic three-cluster state |<?>j) is an eigenstate of the channel Hamiltonian 
Hq + Vj with the eigenvalue JS, i.e., 

\<l>j)=GoVj\$j). (9) 

It is given by the bound state wave function for the pair j times two plane waves 
corresponding to the relative motion of the three free clusters. Since the full wave 
function ([6]) is an eigenstate of the Hamiltonian ([1]), i.e., Yl^i=i '"iWp.n) = GQ^\9p.n), 
we rewrite the amplitude ([8]) as 

{<l>,\T,p\<l>p.n) = (<?j|(l - l.,Go)Go >p,n)- (10) 

Although {$j\{l — VjGo) = due to Eq. one has to keep in mind that t; arising 
from the wave function ^ has a pole for \(l>j) if / = j. Using Eq. ^ one gets the 
relation 

(<?jl(l-t;,Go)t,Go = (11) 
that leads to the three-cluster breakup amplitude 

{'^j\T3p\'l'p,n) = ^{#j|C/f Go tfc GoW;^j,|<^^,„). (12) 

'yki 

Finally, the amplitude for the four-cluster breakup can be found in Ref . [20| . 

In the system of identical particles the number of distinct partitions and AGS 
operator components reduces to two. We choose those partitions to be ((12)3)4 
and (12) (34) and denote them in the following by a = 1 and 2, respectively. The 
amplitudes (|12p after symmetrization become 

{<l'd\Tdc.\'^a,n) = Sdc,{<Pd\[i^ + eP3i)UiGotGoUi» + U2GotGoU2o.Ua,n). (13) 

Thus, Tj^i (2^2) describes the atom-trimer (dimer-dimer) breakup into a dimer and 
two atoms. S^i = v'S and 5^2 = 2 are the corresponding symmetrization factors, 
P34 is the permutation operator of particles 3 and 4, and e = 1 (-1) for bosons 
(fermions). The dimer-atom-atom three-cluster channel state \<P^) is symmetrized 
only within the bound pair (12). There is only one independent Faddeev component 
l^a.n) for each two-cluster channel state. The two-particle transition matrix t is 
derived from the pair (12) potential according to Eq. The symmetrized 3+1 
and 2 + 2 subsystem operators are 



Uc.= PaGo' +PatGoUc 



(14) 
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with Pi = Pi2 P23 + Pi3 P23 and P2 = P13 P24 where P^f, is the permutation 
operator of particles a and b. The symmetrized four-particle transition operators 
obey the symmetrized AGS equations [25], i.e., 

Uii = eP34{GotGor^ + eP34UiGotGoUii + U2GotGoU2i, (15a) 

W21 = (1 + eP34){GotGoy^ + (1 + eP3i)UiGotGoUii, (15b) 

W12 = {GotGor^ + eP34f/iGotGoWi2 + U2GotGoU22, (15c) 

W22 = (1 + eP34)UiGotGoUi2. (15d) 

The employed basis states have to be symmetric (antisymmetric) under exchange 
of two bosons (fermions) in subsystem (12) for the 3 + 1 partition and in (12) and 
(34) for the 2 + 2 partition. 

After partial-wave decomposition for each combination of the total angular 
momentum J and parity 77 the AGS equations (|15(1 become a system of coupled 
three-variable integral equations. In our momentum-space framework these vari- 
ables are the magnitudes of the Jacobi momenta kx, ky and kz that are defined 
in Ref. [20j and describe the relative motion in the 1+1, 2+1, and 3+1 (1+1, 
1+1, and 2+2) subsystems of the 3+1 (2+2) configurations, respectively. Such 
three- variable integral equations have been solved in Refs. [27,28J for various four- 
nucleon reactions with two-cluster initial and final states. However, by using a 
separable two-particle potential the AGS equations (|15p can be reduced to a sys- 
tem of integral equations in two continuous variables, the Jacobi momenta ky and 
kz- We applied this technical simplification in the study of the universal proper- 
ties of the four-boson system that are independent of the short-range interaction 
details |141I171[20] . The same separable potential is used also in the present work. 
The discretization of integrals in the AGS equations using Gaussian quadrature 
rules [29J leads to a system of linear algebraic equations. While the calculational 
technique to a large extent is taken over from Refs. [141I171[20] . there is one im- 
portant difference associated with the open dimer-atom-atom channel. Namely, in 
addition to the subsystem bound state poles of Ua the kernel of the AGS equa- 
tions has integrable singularities arising from t. For ky/2fj,ay + k1/2^a E — 
the two-boson transition matrix in the channel with the dimer quantum numbers 
for the pair (12) has the pole 



iJ + iO — Erf — ky/2flay — ki/2iia 

where e,; < is the dimer energy and ^lay is the reduced mass for ky. The treatment 
of this singularity simplifies considerably if the available four-boson energy is just 
at the dimer-atom-atom threshold, i.e., E = e^. In this case there is only one 
singular point ky = kz = and in all nontrivial momentum integrals (the ones 
without (5-functions) of the type 



L 



f{ky^ kz)kj dk 



^0 - k^/2flay - ki/2fia 



(17) 



this singularity is cancelled by fc| leading to numerically harmless integrals. 
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3 Dimer-atom-atom recombination 

We consider an ultracold gas consisting of a mixture of bosonic atoms and dimers 
with respective particle densities pa and p^. We define the dimer-atom-atom recom- 
bination rate K211 such that the number of the dimer-atom-atom recombination 
events per volume and time is _K'2iiPa/5d/2!. At sufficiently low temperature the 
momenta of the initial particles are much smaller than the final momenta of the 
two resulting clusters, i.e., ky,kz « Pa,n'^ s-nd the kinetic energies of initial par- 
ticles are much smaller than the final two-cluster kinetic or binding energies, i.e., 
k'^/2fj.ay + k^jlpa « Pa,ra/2pa, |ea,n|. Under these conditions one can assume 
the threshold values ky = kz = for the initial three-cluster channel state \<Pj} 
and E = ej_ for the AGS equations. The dimer-atom-atom recombination rate in 
this limit is given by 

KfdlTdal'^a,,.')!'. (18) 

Here we used the time-reversal symmetry in replacing the amplitude for the three- 
cluster recombination into the final two-cluster state by the amplitude for the 
three-cluster breakup of the initial two-cluster state, i.e., 

{^a.n\T^d\^a} = {^d\Tda\^a.n}. (19) 

Furthermore, at fcj, = fcz = only the J' = partial wave contributes to the 
amplitude dH}, i.e., {$d\Tda\^a,n) = {^^V^'^i^diJ = G)\Td^\$a,n{J = 0))- 

The dimer-atom-atom recombination process leads to the final state either 
with two dimers or with an atom and an Efimov trimer. As it is evident from 
Eq. (|18p . K211 can be decomposed into the dimer-dimer and (n -\- 1) atom-trimer 
contributions 

n 

i^2ii =i^2n+ E ^211, (20) 

n'=0 

where n labels the highest excited trimer state. 

We aim to determine the universal behaviour of the ultracold dimer-atom-atom 
recombination rate and its contributions (|20p as functions of the two-boson scat- 
tering length. We therefore build dimensionless quantities a/a^ and K2iim/Tia'^, 
with m being the boson mass. Due to asymptotic discrete scaling symmetry it 
is sufficient to consider the regime 1 < a/ai < a^^-^/ai ^ 22.694 [a with n + 1 
trimer states, provided n is sufficiently large such that non-universal short-range 
corrections become negligible. We found that n > 3, i.e., one has to consider the 
processes involving four or more atom-trimer channels. This is fully consistent 
with previous findings |10II14|[20| for other four-boson reactions. 

We compare the most important contributions to the ultracold dimer-atom- 
atom recombination rate in Fig. [2] All contributions peak at a/a$^ = 1 and a/af^ ~ 
22.694 where the n-th and the (n -|- l)-th trimers, respectively, emerge at the 
atom-dimer threshold. Near these critical regimes the recombination process is 
dominated by the two-cluster channel with the weakest binding, i.e., the n-th atom- 
trimer channel close to a/a^ = 1 and the dimer-dimer channel close to a /an ~ 
22.694 with the respective contributions -^"211 ^t-nd -K^2ii- The n-th atom-trimer 
and dimer-dimer thresholds intersect at a/un ~ 6.789 [IT]- However, it seems that 
there is an interference between these two channels around a/a^ = 10 leading to 
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Fig. 2 (Color online) Dimer-atom-atom recombination rate as a function of the two-boson 
scattering length. The contributions of the final dimer-dimer (dashed-dotted), n-th (solid) and 
(n — l)-th (dashed) atom-trimer channels are compared with the three-atom recombination 
rate (dotted). 

a slight enhancement of K211 and suppression of -R'211 such that K211 > -^"211 
for a/an < 13.5. This behaviour is not very surprising given the strong couphng 
between the n-th atom-trimer and dimer-dimer channels observed in the dimer- 
dimer scattering [17; . The transitions to channels with lower trimers (n' < n — 1) 
are strongly suppressed; we only show the predictions for A'jil^, the others being 
considerably smaller. 

In addition to various K211 contributions we show in Fig. [2] also the ultracold 
three-atom recombination rate K3 defined such that the number of three-atom 
recombination events per volume and time is A'apa/S!. We calculate it numerically 
as 

K3 = 2{27T)\iyPy\{<t>o\{l + Pi)tGoUiMf (21) 

where \ (f>j} and \(j>o) are the atom-dimer and three-atom channel states and py is the 
relative atom-dimer momentum in the final state. Our K3 predictions are in perfect 
agreement with the semi-analytical results of Ref. [3] (after taking into account 
that the respective definitions differ by a factor 3!). We observe that A"2ii is much 
larger than K3, especially near the ends of the shown interval and around the 
minimum. This may have important consequences for the life time and stability 
of trapped ultracold atomic gases with small admixture of dimers: Let's suppose 
that the trapping potential is lower than the kinetic energies of the recombination 
products, i.e., they are able to escape the trap. Neglecting other reactions like the 
four-atom or dimer-dimer-atom recombination that should be suppressed at not 
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too high densities, the atom density in the trap evolves with time according to 

dpa 1 3 2 /r,o\ 

= " 2 - K21lPaPd- (22) 

Thus, due to K211 » even at low dimer densities pj^ « pa the dimer-atom- 
atom recombination may be as important as the three-atom recombination. 

Finally, one may raise the question whether a complicated four-body treatment 
of the dimer-atom-atom recombination is really necessary and to what extent a 
simple three-body model considering the dimer as a point-like particle is reliable. 
Such model was constructed in Ref . [30' . Although it describes well the nth trimer 
binding and the low-energy atom-dimer scattering, it does not support lower atom- 
trimer channels, i.e., for n' < n it yields ifJii = 0. Furthermore, the three-body 
dimer-atom-atom model is obviously inappropriate in the dimer-dimer channel 
since it treats the two dimers asymmetrically. Consistently with two-cluster reac- 
tions from Ref. [30 , the three-body dimer-atom-atom model fails heavily for K211 
as shown in Fig.|3l However, -ffJii in the regime a/a^ < 4 is reproduced quite well, 
with 10% or better accuracy. There the trimer size exceeds the dimer size such that 
the approximations of a point-like dimer and a trimer consisting of weakly bound 
atom and dimer are reasonable. For a/ an > 4 the three-body dimer-atom-atom 
model fails also in predicting -?f2ii- 



4 Summary 

We studied the dimer-atom-atom recombination in the universal four-boson system 
that exhibits the Efimov effect. We used exact AGS equations for the four-particle 
transition operators that were solved in the momentum-space framework. The 
relation between the AGS operators and the dimer-atom-atom recombination am- 
plitude was derived. We calculated the dimer-dimer and atom-trimer contributions 
to the dimer-atom-atom recombination rate in the ultracold limit. We found that 
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usually the channel with the weakest binding dominates the recombination process 
but there are deviations from this behaviour when several two-cluster thresholds 
are close to each other. We also have shown that the dimer-atom-atom recombina- 
tion rate greatly exceeds the three-atom recombination rate leading to important 
consequences for the life time of atomic gases with a small admixture of dimers. 
We also studied a simplified three-body model for the dimer-atom-atom recombi- 
nation process and found it to be inadequate, except for K211 in the regime where 
the trimer can be considered as a weakly bound atom-dimer system. 
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